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Abstract 

The Born effective charge tensors of Barium Titanate have been calculated 
for each of its 4 phases. Large effective charges of Ti and O, also predicted by 
shell model calculations and made plausible by a simplified model, refiect the 
partial covalent character of the chemical bond. A band by band decomposi- 
tion confirms that orbital hybridization is not restricted to Ti and O atoms 
but also involves Ba which appears more covalent than generally assumed. 
Our calculations reveal a strong dependence of the effective charges on the 
atomic positions contrasting with a relative insensitivity on isotropic volume 
changes. 

PACS numbers: 77.84.Dy, 77.22.Ej, 71.25.Tn 
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Barium Titanate (BaTiOs) is a well known ferroelectric compound |T|]. Its structure, 
cubic perovskite at high temperature, becomes tetragonal around 120°C, orthorhombic at 
about 0°C and rhombohedral near -70°C. Although it is probably one of the most studied 
ferroelectric crystal ||1||, the mechanism of its phase transitions is still subject of controversy 
0. The character of the chemical bonding also remains questionable 

The Born effective charge Z* is a fundamental quantity for the study of lattice dynamics, 
controlling the long range Coulomb part of the force constants. An estimation of Z* for 
three ABO3 compounds was first proposed by Axe on the basis of a phenomenological 
theory. However, the crude hypothesis inherent in his procedure only allowed an approximate 
estimation of these charges and restricted his investigations to the cubic phase. New advances 
in "ab initio" techniques now enable one to determine Z* theoretically using perturbation 
theory ||^ or finite difference of polarisation . Accurate values have been recently reported 
by Resta et al. for KNbOs 0, by ourselves for BaTiOa and by Zhong et al. for a 
whole set of ABO3 compounds. These studies underline the surprisingly large values already 
observed by Axe and generally referred to as "anomalous" charges. For BaTiOa, no attempt 
was performed to determine how these giant charges are affected by structural details, even 
if this information could reveal essential for an accurate investigation of lattice dynamics. 
We will show that Z* are relatively insensitive to isotropic volume changes but are strongly 
affected by changes of positions associated to the phase transitions. 

Moreover, if obtaining the values of Z* is already an important result, to understand why 
some charges are so large also constitutes a fondamental problem. Until now, anomalous 
charges in ABO3 materials were explained qualitatively in the framework of a bond orbital 
model, following Harrison [0. Recently, Posternak et al. [|lT| elegantly emphazised the role 



of covalency between Nb and O in KNbOs from the analysis of giant effective charges. We 
show here that, in a more general way, a band by band decomposition of Z* is a sensitive 
tool to investigate the role of covalency and ionicity without any preliminary hypothesis on 
the orbitals that interact. Our study helps to clarify the debate on the real nature of the 
chemical bonding in BaTiOs and brings out the role of the Ba atom. 



We work in the framework of the Density Functional Formahsm within the local density 
approximation . The Born effective charge tensor Z* of atom k can be linked either 
to the change of polarisation induced by the periodic displacement r^^a, or to the force 
induced on atom k by an electric field 8^: Z*^ = = ^fg'" = —j^MP- — . So, 

it also appears as a mixed second derivative of the total energy E per unit cell volume V, 
which can be evaluated using the variational approach to Density Functional Perturbation 
Theory We use the exchange-correlation energy functional and the norm- conserving 

pseudopotentials of ref. pf. The Brillouin Zone was sampled with a 6x6x6 Appoint set. 
The plane-wave basis set was determined by a 35 hartrees energy cutoff which guarantees a 
convergence better than 0.5% on Z*. 

First we investigate the Born effective charges in the cubic structure. Ab-initio values, 
obtained at the experimental (aceii=4A) and theoretically optimized (acez/=3.94A) lattice 
parameter, as well as those corresponding to a compressed cubic cell (ace//=3.67A) are re- 
ported in Table |. In comparison, we also computed the effective charges from the shell 
model parameters proposed by Khatib et al. ||13|| . 



The charges of Ba and Ti, isotropic owing to symmetry, are equal respectively (at the 
theoretical volume) to +2.77 and +7.25. For O, the values of (-5.71) and Zq^ (- 
2.15) refer respectively to a displacement of the oxygen ion along the Ti-0 direction or 
perpendicular to it. Our results are in good agreement with those obtained by Zhong et al. 
(also at the theoretical volume) from finite difference of polarisation and globally reproduce 
the values deduced by Axe ||^. The accordance with the shell model is only qualitative 
and illustrates the limited precision obtained within such an empirical description. We note 
however that this shell model was not designed to reproduce the dielectric properties of 
BaTiOa. The large value of Z^^ (7.51) proves that it implicitely includes covalency effects 
responsible for anomalous effective charges, as described below. 

The charges on Ti and Oy are surprisingly large in the sense that they reach about 
twice the value they would have in a pure ionic picture: they reveal the presence of a large 
dynamic contribution superimposed to the static charge. As the latter quantity is ill defined. 



determining the dynamic contribution is an ambiguous task. So, in the following, we prefer 
to speak in terms of anomalous contributions that we define as the additional charge with 
respect to the well known, ionic value (cfr. first column of Table 1). 

The physical possibility of obtaining anomalous charges can be understood within a 
very simplified model. Let us just consider a diatomic molecule XY with an interatomic 
distance u and a dipole moment P{u). The dipole moment enables us to define a static 
charge Z{u) = and a dynamic charge Z*{u) = ^^^^ also equal to Z{u) + u^^^. As 
the distance between X and Y is modified from to some u (the distance corresponding 
to a complete transfer of electrons from X to Y), the dipole moment evolves continuously 
from P{0) = (since there is no dipole for that case) to P(u). Interestingly, Jq Z*{u)du = 
[P{u) - P(0)] = uZ(u). So, i J^Z*{u)du = Z{u): the mean value of Z*{u) from to U is 
equal to Z(u). This result guarantees that, if Z{u) changes with u, Z*{u) has to be greater 
than Z{u) for some u between [0,ll]. The difference between Z{u) and Z*{u) can be large if 
Z{u) changes rapidly with u. 

In BaTiOa, the approximate reciprocity between 0_l (-3.71) and Ti (+3.25) anomalous 
contributions suggests that they correspond to a global transfer of charge from O to Ti 
when the Ti-0 distance shortens. In the framework of the bond orbital model proposed 
by Harrison |jlO|, the charge redistribution is attributed to a change in the hopping integral 



produced by dynamic modification of orbital hybridizations. Matching this bond orbital 
model to the real material ask for the identification of the relevant orbitals. 

Hybridization between O 2p and the unoccupied metal d orbitals is well known and was 
already pointed out from experiments [|1^, tight-binding models [|10| and LCAO calculations 



||15|| . It was highligthed more recently by Cohen from first-principles as an essential 
feature of ABO3 ferroelectric compounds. In this context, it seemed realistic, following 
Harrison, to focus on O 2p-Ti 3 d hybridization changes to explain intuitively large anomalous 
contributions |^. Posternak et al. went beyond this credible assumption. They showed 
for KNbOs that the anomalous contributions disappear when the interaction between the 
O 2p and Nb Ad orbitals is artificially suppressed. 



Nothing proves that the conclusions would be the same for BaTiOa neither that the 
hybridizations are limited to this kind of p-d interactions. Theoretical investigations, con- 



firmed by experimental results suggest that Ba also plays a major role in forming the 



valence band structure. Now, we propose a more detailed investigation, based on a band by 



band analysis of Z* 



Contrariwise to the total effective charge tensor, the band by band decomposition de- 
pends not only on the Hilbert space of occupied valence states but also on the particular 
valence eigenfunctions. In order to identify the contribution of each band we worked in the 
"diagonal gauge" [|l^] in which the matrix of the first-order eigenenergies with respect to 



atomic displacements is diagonal. This is equivalent to associate Wannier functions with 



each separate set of bands for calculating the polarization ||6|JT8[ . 

For a reference configuration in which the contribution to the charge of a given atom 
is -2 for each band associated to its own orbitals and for the other bands, each Wannier 



function is centered on a given atom. Within the recent theory of the polarization PJT8|], the 
anomalous charge associated to a particular band for a given ion (defined as the additional 
part with respect to our reference value) reflects how the center of the Wannier function 
of this band moves with respect to the ion. In a purely ionic crystal, each band would be 
composed of non hybridized orbitals. As covalency develops, different orbitals mix adding 
anomalous contributions from several bands. 

Results for the theoretical cubic structure are reported in Table 0. The first line brings 
together the charge of the nucleus and the core electrons included in the pseudopotential. 
The other contributions come from the valence electron levels. Sets of bands were identified 
by the name of the main atomic orbital which generated this energy level in the solid (Fig. p. 
Their dominant character was confirmed by partial density plots. As a check of the coherency 
of our decomposition, the sum of the contributions of each fully occupied band, obtained 
independently, is -2 (within ±0.01, the accuracy of our calculation). 

As expected, the main Z^^ anomalous charge is localized on the O 2p bands (+2.86). It 
can be understood by an hybridization between O 2p and Ti 3 (i orbitals. Interestingly, there 



are also smaller but non-negligible anomalous charges from the Ti 3p (-0.22), O 2s (+0.23) 
and Ba 5p (+0.36) bands. The different positive contributions correspond to a displacement 
of the center of the Wannier function of the O and Ba bands in the direction of Ti when 
this atom moves. The Ti3s contribution is close to -2. This result a posteriori justifies the 
inclusion of deeper levels as part of ionic pseudopotentials. 

For Z^^, the decomposition is more surprising: the small global anomalous effective 
charge (+0.77) that could be typical of a more ionic character appears partially fortuitous. 
The anomalous charges of O 2s (+0.73) and O 2p (+1.50) bands are not small. Nevertheless, 
they are partially compensated by Ba 5p (-1.38) and Ba 5s (-0.11) anomalous contributions. 
This result gives a tangible proof of the hybridization of Ba 5p orbitals already suggested 



from the experiment [14 



Concerning oxygen, even if Oy and are defined respectively for a displacement of O 
in the Ti and Ba direction, it seems only qualitative to associate with and Zq^ with 
Z^^. The O 2p anomalous contributions to Ti and Oy do not compensate. Moreover, O 2p 
contribution to Z'^^ not only comes from but has equivalent contributions coming from 
0||. 

Within this analysis, several bands appear as a complex mixing of orbitals coming from 
the different ions. In this context, and contrasting with the conclusion of Posternack et al. for 
KNbOa, a correct understanding of the Born effective charge goes here beyond the simple 
model of Harrison. Our result clarifies the mixed ionic-covalent character of BaTiOs: it 
clearly establishes that the covalent character is not restricted to the Ti-0 bond as generally 
assumed. 

Until now, calculations of Z* essentially focused on the cubic phase P,^. On the basis of 
an investigation of these charges in the experimental tetragonal structure of KNbOs [0 and 
PbTiOa , it was argued that they are quite insensitive to structural details. This result is 
surprising if we remember that anomalous contributions to Z* are closely connected to orbital 
hybridizations, these in turn, being strongly affected by phase transitions In addition, the 
theoretical overestimation of the spontaneous polarization for the rhombohedral structure 



of BaTiOs also suggests a reduction of Z* in this phase. 

We computed the Born effective charge tensors for the 3 ferroelectric phases at the 



experimental unit cell parameters |J19|], with relaxed atomic positions |2y]. In this paper, 
we only comment on the eigenvalues of these tensors (Tables pn] - |IV|) that already allow a 
pertinent comparison with the cubic phase. The Z33 eigenvalues of Ba and Ti correspond 
to an eigenvector aligned along the ferroelectric axis. In the case of O, the eigenvector 
associated to the highest eigenvalue approximately point in the Ti-0 direction: we identify 
this highest contribution as Oy while the others are refer as 0_l, by analogy with the cubic 
phase. 

Although the charges of Ba and remain globally the same in the 4 phases, for Ti and 
0||, stronger modifications are observed. Changing the Ti position by 0.076A (2% of the unit 
cell length) when going from the cubic to the rhombohedral phase, reduces the anomalous 
part of by more than 50% along the ferroelectric axis (Table p|) . The amplitude of 
and Zq in the direction of the shortest Ti-0 bond length dmin of each phase is plotted 
in Fig. with respect to the interatomic distance dmin- For the different phases at zero 
pressure, the anomalous parts decrease with dmm- The comparison with a compressed cubic 
phase at 3.67 A in which the Ti-0 distance is the same that the shortest Ti-0 bond length 
in the tetragonal structure shows nevertheless that the evolution of Z^^ cannot be explained 
in terms of the Ti-0 distance only but is critically affected by the anisotropy of the Ti 
environnement . 

This is clear from the band by band decompositions of in Table ||. While in the 
cubic structure at 3.67 A every Ti-0 distance is equivalent to the others, in the tetragonal 
phase, along the ferroelectric axis, a short Ti-0 bond length is followed by a larger one which 
breaks the Ti-0 chain in this direction and inhibits the giant current associated to the large 
effective charges. This appears at the level of the O 2p bands (-1-1.48 instead of -f 2.86) while 
the other contributions remains equivalent to those of the cubic phase at 3.94 A. Analysis 
from the cubic structure at 3.67 A reveals that the O 2p contribution is not significantly 
affected by hydrostatic pressure; on the other hand, the anomalous parts of the Ba 5p, Ba 



5s and Ti 3p bands are modified by about 50 % due to the compression. 

In this work, we were able to compute the Born effective charges of BaTiOa in its 4 
phases. Effective charges are a sensitive tool for analysing dynamic changes of orbitals 
hybridizations, especially if a band by band decomposition is performed. In our description 
Ba appears more covalent than generally assumed. The charges of Ti and O are strongly 
affected by atomic displacements but quite insensitive to hydrostatic pressure. 
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from common projects between IBM-Belgium, UCL-PCPM and FUNDP. 
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FIGURES 

FIG. 1. Electronic band structure of cubic BaTiOs (aceii=3.94A). 

FIG. 2. Born effective charge of Ti (opened symbols) and O (filled symbols) in the direction of 
the shortest Ti-0 bond length dmin with respect to this interatomic distance for the cubic (square) , 
tetragonal (lozenge), orthorhombic (circle) and rhombohedral (triangle) phases. 
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TABLES 









TABLE I. Born effective charj 


2;es for cubic 


BaTiOs. 






Ions 


Axe [| 


1 Shell Model 


Zhong et al. | 




!=4.00A acei«=3.94A 




7* 
^Ba 


+2 


2.9 


1.63 


2.75 




2.74 


2.77 


2.95 




+4 


6.7 


7.51 


7.16 




7.29 


7.25 


7.23 




-2 


-2.4 


-2.71 


-2.11 




-2.13 


-2.15 


-2.28 


z* 




4 8 


-3.72 


-5.69 




-5.75 


-5.71 










TABLE IL Band by band decomposition 


of Z* (see text). 










cubic(3.94A) 






cubic(3.67A) 


tetragonal 


Band 




^Ba 














core 




iU.UU 


6.00 


6.00 


12.00 




12.00 


1 o nn 
iZ.OO 


il OS 




U.Ui 


0.00 


0.02 


-2.03 




-2.07 


o nc^ 
-2.05 


Ti 3p 




0.02 


-0.02 


0.21 


-6.22 




-6.43 


-6.26 


Ba 5s 




-2.11 


0.02 


0.01 


0.05 




0.09 


0.05 


2s 




0.73 


-2.23 


-2.51 


0.23 




0.27 


0.25 


Ba 5p 




-7.38 


0.58 


-0.13 


0.36 




0.64 


0.34 


2p 




1.50 


-6.50 


-9.31 


2.86 




2.73 


1.48 


Total 




2.77 


-2.15 


-5.71 


7.25 




7.23 


5.81 
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TABLE III. Eigenvalues of the Born effective charge tensors of Ba and Ti at the experimental 

volumes. 



phase 














^Sa,22 


y * 

■^Sa,33 


cubic 


7.29 




7.29 


7.29 




2.74 




2.74 


2.74 


tetragonal 


6.94 




6.94 


5.81 




2.72 




2.72 


2.82 


orthorhombic 


6.80 




6.43 


5.59 




2.72 




2.81 


2.77 


rhombohedral 


6.52 




6.52 


5.58 




2.79 




2.79 


2.74 


TABLE IV. Ei| 


^envalues of the Born effective charj 


;e tensors of at the experimental volumes. 



















0|| 




phase 


rr* 

^oi,ii 


^02,11 


^03,ll 


ry* 

^01,22 


^02,22 


^03,22 


y* 

^01,33 


y* 

^02,33 


y* 

^03,33 


cubic 


-2.13 


-2.13 


-2.13 


-2.13 


-2.13 


-2.13 


-5.75 


-5.75 


-5.75 


tetragonal 


-1.99 


-1.95 


-1.95 


-1.99 


-2.14 


-2.14 


-4.73 


-5.53 


-5.53 


orthorhombic 


-1.91 


-1.91 


-1.97 


-2.04 


-2.04 


-2.01 


-4.89 


-4.89 


-5.44 


rhombohedral 


-1.97 


-1.97 


-1.97 


-1.97 


-1.97 


-1.97 


-5.03 


-5.03 


-5.03 
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Ghosez et al., Fig. 1 
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